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The Random Transverse Field Ising Chain is the simplest disordered model presenting a quan- 
tum phase transition at T = 0. We compare analytically its finite-size scaling properties in two 
different ensembles for the disorder (i) the canonical ensemble, where the disorder variables are 
independent (ii) the microcanonical ensemble, where there exists a global constraint on the disorder 
variables. The observables under study are the surface magnetization, the correlation of the two 
surface magnetizations, the gap and the end-to-end spin-spin correlation C{L) for a chain of length 
' L. At criticality, each observable decays typically as e"™^ in both ensembles, but the probability 

' distributions of the rescaled variable w are different in the two ensembles, in particular in their 

, asymptotic behaviors. As a consequence, the dependence in L of averaged observables differ in the 

' two ensembles. For instance, the correlation C{L) decays algebraically as 1/L in the canonical en- 

^ ^. semble, but sub-exponentially as e"'^^ in the microcanonical ensemble. Off criticality, probability 

■ distributions of rescaled variables are governed by the critical exponent = 2 in both ensembles, 

' but the following observables are governed by the exponent D = 1 in the microcanonical ensemble, 

instead of the exponent v = 2 in the canonical ensemble (a) in the disordered phase : the averaged 
surface magnetization, the averaged correlation of the two surface magnetizations and the averaged 
end-to-end spin-spin correlation (b) in the ordered phase : the averaged gap. In conclusion, the 
measure of the rare events that dominate various averaged observables can be very sensitive to the 
microcanonical constraint. 
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I. INTRODUCTION 
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CO , A. Microcanonical ensemble versus canonical ensemble for a finite disordered sample 

o 

i In the study of disordered systems, it is usual to consider that the random variables defining the disorder in a given 
■ sample are independent : following [3, Ql j this procedure will be called here the "canonical ensemble" (this procedure 
^ , is also called the "grand-canonical ensemble" in However, it has been argued in that it is much more 

' interesting in some cases to consider the so called "microcanonical ensemble" as in 0, ( this procedure is called 
^ ' the "canonical ensemble" in H, ^ |M IS) • in the microcanonical ensemble, there exists a global constraint on the 
O random variables defining the disorder in a given sample of N sites. The first important example concerns a pure 
system with a fraction p € [0, 1] of impurities fsj, i.e. the disorder is characterized by a binary distribution : in the 
canonical ensemble obtained by putting independently on each site an impurity with probability p, the total number 
of impurities presents fluctuations of order ^/N around its mean value pN, whereas in the microcanonical ensemble, 
the total number of impurities is fixed to be exactly pN and presents no fluctuations at all, i.e. the remaining disorder 
only concerns the positions of the impurities. It is thus clear that the microcanonical ensemble is much less disordered 
that the canonical ensemble, and indeed, it may yield a spectacular noise reduction in numerical studies 3, 4]. Another 
example where the ensemble dependence has been under study recently is the Random Transverse Field Ising Chain 
P, , deflned by the Hamiltonian 

i/ = -^J.afaf+i-^/i.ar (1) 

i i 

where the couplings Ji > and the fields hi > are random variables (the signs can be fixed via a gauge trans- 
formation). This model is of great interest because it is the simplest disordered model presenting a quantum phase 
transition at T = : the critical point situated at 

[\llh]av = [In J]a« (2) 

separates a ferromagnetic phase ([In/iJ^^ < [In J]a„) from a disordered phase ([In/ija^, > [In J]q„). This model can be 
studied in great details via a disorder-dependent real space RG analysis Q , which agrees with previous exact results 
whenever they can be compared. In numerical studies of the quantum critical point on systems of linear size 
L, the usual canonical ensemble has been compared to the microcanonical ensemble defined by the global constraint 
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^(lnJ,-ln/iO =0 (3) 

i=l 

(with one fixed boundary condition so that there are exactly the same number L of random bonds and random fields) . 
In the canonical ensemble, the quantity on the left in equation Q presents fluctuations of order \/L around its mean 
value zero. The samples of the microcanonical ensemble are thus 'closer to criticality' in some sense than the samples 
of the canonical ensemble. So in both examples, the microcanonical constraint appears as an opportunity to obtain 
numerical results with much less fluctuations. 

Now the question is : is that interesting from the physics point of view? On one hand, it has been strongly argued in 
0,3 that the microcanonical ensemble should be preferred to the canonical ensemble, because the latter introduces an 
"extra noise" that may hid the "intrinsic" properties of the system. On the other hand, it seems that the fluctuations 
of order -s/I for the sum of I random variables are precisely an essential property of disordered systems, and indeed 
these fluctuations appear in the Harris argument [13] on the relevance of disorder near critical points, in the Imry-Ma 
argument for random field Ising models, and in the Chayes et al. theorem for the bound i' > 2/d. Moreover, 
if one divides the system of size L into two halts of size L/2, each half will present fluctuations of order y/T in both 
ensembles : in the canonical ensemble, these two halts are independent, whereas in the microcanonical ensemble, 
the two halts are completely correlated, i.e. they have exactly opposite fluctuations. From this point of view, the 
microcanonical constraint can thus appear to be quite artificial or even biased. 

Of course, it seems a priori natural to expect that these two ensembles should be equivalent in the thermodynamic 
limit, as was shown in 6] for the case of a random classical ferromagnet. However, it is clear on the two examples 
above that their flnite-size properties can be very different. But since the finite-size scaling theory of phase transitions 
relates the thermodynamic exponents to flnite-size effects obtained in numerical simulations, the discussion about 
these two ensembles actually leads to the general problem of the finite-size scaling theory for disordered systems 



B. Relation with the finite-size scaling theory for disordered systems 

To generalize the ideas of the finite-size scaling theory developed for pure systems to the case of disordered sys- 
tems, the central question is |^ |^ |^ : when the numerical simulations provide values for some observable, like the 
susceptibility x(i,L), for various samples (i) and various sizes L, what is the best procedure to analyze these data? 
The usual procedure consists in averaging the observable Xm{L) —< x(i,L) >i over the samples (i) at fixed size L 
and in analyzing the dependence in L of the mean values Xm{L) as in the finite size scaling theory of pure systems. 
However, as explained in details in (31,(3, Is'l with illuminating numerical examples, the correct procedure to analyze the 
data consists in a finite-size scaling analysis for each given sample. In particular, one has to defined a pseudo-critical 
temperature Tc(i, L) for each sample, either by the maximum of the susceptibility 0,(3 or by a distance minimizing 
procedure P. Indeed, these studies have shown that the fluctuations of xihL) over the samples are mainly due to 
fluctuations in Tc{i, L), and that the use of the variable (T — Tc{i,L)) in the flnite-size scaling analysis yields much 
more accurate results (see for instance Fig. 2 and Fig. 3 in Q). Within this point of view, the advantage of the mi- 
crocanonical ensemble with respect to the canonical ensemble is justified if the microcanonical constraint determines 
Tc{i, L) or at least reduces drastically its fluctuations over the samples, as it is the case for binary distributions H,^. 
For quantum phase transitions, the pseudo-critical temperature T(.{i^ L) of classical systems has to be replaced by the 
pseudo critical point 4]. In conclusion, the best justiflcation for the microcanonical ensemble seems to be in systems 
where the microcanonical constraint actually corresponds to a constraint on the pseudo-critical temperature like in 
critical phenomena with a binary distribution (3, or to a constraint on the pseudo critical point for quantum phase 
transitions as in the Random Transverse Field Ising model. 



C. Ensemble dependence in the RTFIC 

The ensemble dependence in the RTFIC has already been studied in 0, |2| . The flrst considered observable has 
been the surface magnetization 7715 : in particular, its averaged value over the samples was found to present different 
scaling in L at criticality in the two ensembles 0, 13I ■ it decays as 1/ Vl in the canonical ensemble and as 1 /L in the 
microcanonical sample However, the interpretation that should be given to this difference is still under debate 

Q Q ' in Q it was interpreted as the presence of the two different correlation exponents v = 2 and = 1 j whereas 
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in [3 it was interpreted as the same exponent v = 2 in both ensembles, with an accidental vanishing of an amplitude 
in the microcanonical ensemble. 

In 0, other observables have been numerically compared at criticality in the two ensembles, from the point of view 
of probability distributions as well as averaged values. The considered observables are the middle-point magnetization, 
the end-to-end spin-spin correlation 

C{L) ^< a{al > (4) 

the gap, i.e. the energy difference between the two lowest levels 

A{L) = Ei-Eo (5) 

and also the correlation (mf to£) between the two surface magnetizations where each surface magnetization corre- 
sponds to the magnetization at one end when the spin at the other end is fixed 

ml =< o\ > \ai=i 

ml =< al > Uj-=i (6) 

We refer to for a detailed discussion on the numerical results for each observable. 

The aim of this paper is to compare analytically the behavior of these various observables in both ensembles, at 
criticality as well as off-criticality, with a particular attention to the averaged values, that are governed by rare events. 

The paper is organized as follows. In Section^ we recall the relation between the surface magnetization and Kesten 
random variables, and introduce the important notations. In Section IIIII we show that the limit distributions for 
the surface magnetization obtained in can be understood as resulting from a saddle point method in each sample. 
In Section llVI and Ivl we compute the exact distribution and the exact averaged value of the surface magnetization, 
in the canonical and microcanonical ensembles respectively, via path-integral methods. In Section IVII we apply the 
saddle-point method to compute the limit distributions for the correlation of the two surface magnetizations in the 
two ensembles. In Section fVIII we compute via the real-space renormalization approach the probability distribution of 
the gap and of the spin-spin end-to-end correlation in the microcanonical ensemble and compare with the previously 
known results for the canonical ensemble Finally, the Section IVlIII contains the conclusions and the appendixes 
some technical details. 



II. SURFACE MAGNETIZATION AND KESTEN RANDOM VARIABLES 



In this Section, we recall that the surface magnetization is a Kesten random variable, and introduce the important 
notations. 



A. Surface magnetization in each given sample 

Remarkably, the surface magnetization © in any given sample of size (L -I- 1) has the following exact expression 

m 

mf = [1 + Zl]-'/^ (7) 

where 

1=1 j=i 



has a specific structure of a sum of products of random variables : we will call a Kesten random variable |l4| . 
Thi s type o f random variables also appears in the context of random walks in random media for various observables 
Ra, uQ, IttL , as well as in the random field Ising chain via the formulation with 2x2 random transfer matrices 
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B. Exponent fj, for discrete Kesten random variables 

It is actually convenient to rewrite the Kesten random variable lO with a varying left boundary 

b i 

i=a j=l 

where yj = [hi/Ji)^ arc the independent random variables. The fundamental property of Z{a,b) is the recurrence 
equation 

Z{a,b)=ya[l + Z{a+l,b)] (10) 

where the random coefficient ya appears multiplicatively : Z{a,b) is thus a multiplicative stochastic process. 

One of the main outcome of the studies on the random variable Z{a, b) is that, in the limit of infinite length 
L = b — a —I- oo, there exists a limit distribution Poo(Z) if 

[lny]a. <0 (11) 

Moreover, the limit distribution then presents the algebraic tail [13. IT^ I20II2H 

where the exponent /i is defined as the positive root /i > of the equation 

nv = 1 (13) 

In the field of random walks in random media, this exponent fi is known to govern the anomalous diffusion behavior 
a; in the domain < /i < 1 0,0, 113, El- In the context of the RTFIC, the exponent ^ defined by 1^311 has 
for analog the RG-invariant exponent — (2A) |2J| defined by [{J/h)'^'^]av — 1, and z = 1/(2A) can be interpreted as a 
dynamical exponent. In the vicinity of the critical point A = 0, one may perform a series expansion in A, to obtain 
the solution as A = 6 + 0{S'^), where 

j; _ [lnf]a« _ [In h]av -[Iri J]av ^^^^ 



[(In - ([In ^^a^ln h) + variln J) 

is the parameter introduced in to measure the deviation from criticality in the real-space renormalization approach. 
The expression (|14ll . which comes from an expansion in the two first cumulants of the variable ln(J/ft,), is thus exact 
if the variable ln( J//i) is Gaussian, and it is only an approximation near the critical point for all other distributions 
which present higher order cumulants. 



C. Continuous version of Kesten random variables 

The continuous version of the Kesten random variable © is the exponential functional 

Z[a,b] = / dxe- ^^'^y^^y'^ (15) 

J a 

where {F{x)} is the random process corresponding to the random variables (— Ini/i) = —2 h\{hi/ Ji) in the continuous 
limit. The analog of the recurrence equation H10|l is the stochastic differential equation 

daZ[a,b]= F{a)Z[a,b]-l (16) 

where the random process F{x) appears multiplicatively, in contrast with usual Langevin equations where the noise 
appears additively. In the limit of infinite length L = b — a —> 00, the condition (|11() to have a limit distribution 
Pao{Z) becomes a condition on the mean value of the process F{x) that should be strictly positive 

Fo = [F{x)U > (17) 
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The exponent /i (|12|l is now determined as the root of the equation H13() 

[e-^'-fo''y^^y\^ = i (18) 

for arbitrary x as long as the process F{x) has no correlation. 

It is interesting to note that the exponential functional ifT^ actually determines the stationary flux [2^ 
that exists in a given Sinai sample [0,L] between two fixed concentration cq and cat = (i.e. particles are injected 
via a reservoir at x = and are removed when they arrive at the other boundary x = N) : it is simply given by the 
inverse of the variable Zl = Z[Q, L] 

Jl^^ (19) 

In some sense, it is the simplest physical observable in the Sinai diffusion, as the surface magnetization is the simplest 
order parameter in the RTFIC : both can be expressed in a simple way in terms of the Kesten random variable of the 
sample. 



D. Case of a Brownian process 



The simplest process for F[x) is of course the case where F{x) is a biased Brownian motion 

< Fix) > =Fo (20) 
< F{x)F{x') > -F^ = 2ad{x - x') (21) 



In this case, the exponent fi solution of l(TF)l reads |2a.l25| 

A* = - (22) 

(T 

Here fi exactly coincides with {—26) even away from criticality , as a consequence of the Gaussian distribution. The 
Brownian process actually corresponds to the fixed point of the real-space renormalization approach Q and can thus 
be used to study the universal properties near the critical point. 



III. LIMIT DISTRIBUTIONS FOR THE SURFACE MAGNETIZATION IN THE LARGE L LIMIT 

In this Section, we discuss the universal limit distributions of the random variable (Inms(L)) in the large L limit 
near the critical point, which have been obtained in ^2.) via an approximation of the recurrence relation H1Q(I for InZ 
by an effective biased random walk with a reflecting boundary at the origin. Here we show that these analytic results 
|3 can actually be understood as resulting from a saddle point method in each sample. This point of view gives a 
clearer control on the validity of the approximation, and yields a better understanding of the limitations of the results 
obtained for the averaged surface magnetization that is governed by rare events. Moreover, this approach can be then 
generalized to obtain results for more complex quantities such as the correlation of the two surface magnetization as 
discussed in Section IVll 



A. Saddle point method in each sample 

In the case of a Brownian process (|21ll , the continuous version of the Kesten process (|15|) can be written as 

Zl= f dxe-^(^) (23) 

where the potential U{x) = dyF{y) is a random walk of bias Fq. At the critical point Fq — 0, this potential thus 
presents fluctuations of order ^/L on the interval [0, L]. As a consequence, it seems natural to evaluate this integral 
in the large L limit by the saddle-point method 

Zl ^ e^^Zyalley (24) 
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where {—El) < is defined as the minimum reached by the process U{x) on the interval [0, L], and where Zyaiiey 
represents the partition function of an infinitely deep Brownian valley, whose probability distribution was studied in 
[28l |. At the critical point, the scaling El ^ \IL shows that there will exists a limit distribution for (lnZL)/\/Z ~ 
lnm|/\/i, and that this limiting distribution is given by the limit distribution of El/VL- In the following, we study 
the distribution of the minimum El at criticality and off criticality. 

To study the probability distribution in the two ensembles, we thus need the following standard result : the 
probability for the biased random walk to go from f7o to U during L in the presence of an absorbing boundary at 
U = Ua reads from the methods of images 



G\u::.oo,iU,Lm = -_ e ^ _e--^-o-»^e (25) 



The joint probability of the end-point U and of the minimum value Umin = —E < when starting at C/q = then 
reads 



(Fo) 



Ua=-E 



(U + 2E) Pf) n (U + 2E~Ff,L)^ 

= e{E > 0)eiU > -E) ;^l^^l, e-^^e ^ (26) 

B. Distribution of El in the canonical ensemble 

In the canonical ensemble, we need the partial law of E when integrating over the end-point U in H26fl 

Pi^''\E)\ = [ dUP[^"\u,E) 

\ cano J 

^eiE>0)e-^'= dl-^-^^e-^ (27) 



The final result reads after an integration by part and with the notation Fq = (j/i (|22() 
Vr^iE) - e{E > 0) _^e-^^^ + Me-''^ [^^ ^e" 



(28) 



which exactly coincides with the formula (10) in with the correspondence /i = —25 and / — aL. Let us first consider 
the limit L — > oo with /Lt fixed : 

For /i > 0, there exists a limit distribution for E which is simply exponential 

rt^o°HE) ^ 9iE > 0)^ie-^^^ (29) 



This distribution actually corresponds to the distribution of barriers against the bias [22l, l29l | in a Sinai potential. 
For /i < 0, the appropriate rescaled variable reads 

8 . (^^^ (30) 

For /X < 0, the rescaled variable £ is distributed with the full Gaussian distribution 

Fi^<°)(f) = ^e"^^ (31) 

whereas for fi — 0, the rescaled variable is distributed with the half Gaussian distribution 

P^!ii:'Hn = 9i£>0)^e''' (32) 
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C. Distribution of El in the microcanonical ensemble 

In the presence of the microcanonical constraint (|65() . which fixes the end-point value of the potential U{L) = FqL, 
we need the distribution of E for the trajectories having exactly U = FqL, i.e. we obtain from the joint distribution 
(|25|l with the notation Fq — fia 

Jrmcro j dE P^-^"' {U = Q , E) V ^T-^ / 

which exactly coincides with the formula (11) in with the correspondence ^ = —25 and I = aL. We now consider 
the limit L — > oo with /i fixed, to compare with the results of the canonical ensemble. 
For /i > 0, there exists a limit distribution that exactly coincides with H29|) 

vi^'.^^liE)^^ e{E > 0)t,e''^^ (34) 



At the critical point ^ = 0, the appropriate rescaled variable is again £ = E j^J AaL, but the corresponding 
distribution 

Pi^rliS) = 0{£ > 0)8£e-^^' (35) 

does not coincide with the canonical result (|32|l . In particular, it vanishes linearly as £ — > 0, in contrast with the 
canonical result which presents a finite density at the origin. 

For /i < 0, the appropriate rescaled variable is not £ (|30|l in contrast with the canonical case, but it is 

v = E-\^i\(7L (36) 

which is asymptotically distributed with the exponential distribution 

^ilc.o^M ^d{v>0)^e-^^ (37) 

which is exactly the distribution of barriers against the bias (|29I34|I : this shows that the statistics for the barrier 
E = U{0) — U{xo) is actually completely determined in this case by the statistics of the barrier v = U{L) ~ U{xo) 
distributed with (|37|l . since they are completely correlated via the microcanonical constraint U{L) — U{0) = v ~ E ~ 
FoL (EH). 

D. Discussion 

We first have to discuss the validity of the saddle-point method. At the critical point /i = 0, where the random 
variable E^ scales as \/L, and in the phase /i < 0, where the random variable El is of order L, the saddle-point 
method is thus well justified at large L, at least for typical samples. On the contrary, in the phase fi > 0, where the 
random variable E^ remains finite, the saddle-point method will be a good approximation only if E is large, i.e. in 
the limit /i — s- near the critical point. 

On the other hand, since whenever the saddle-point analysis is valid, the variable El is large, we may safely 
approximate the surface magnetization by 

Inml = \n{l + Zl)-^/^ =i -(1/2) InZ^ - (38) 

Zl>>1 

and thus the results derived above for the distribution of El indeed represents the probability distribution of the log 
of the surface magnetization. In particular, in the vicinity of the critical point, if one introduces the scaling variables 

HQ 



VctL V (jL 

7 = ^iV^ (40) 
we may write the following finite-size scaling forms for the probability distribution of (— Inmf ) 

PLi-lnml) = -^Q (w = ^^i^;7 = (41) 
vaL \ VaL J 
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where the scaling functions respectively read for the two ensembles (|28I33|I 

f +00 



Qcano{w;i) = 0{W > 0) 



4 +76 / -^e 



(42) 



Qm^croiw■,J) = > O)0(w > -7) (2u; + 7) e"^'"""' (43) 

As emphasized in the critical exponents are the same in the two ensembles (the scaling variables (^,7) are the 
same), but the finite-size scaling functions are different in the two ensembles. 

We now turn to the question of the mean surface magnetization. As stressed in 0, Q , since the typical surface 
magnetization decays in the critical region as 

m^=e-'^'^' (44) 

where w is a random variable of order one, the mean surface magnetization will be governed by the rare samples that 
presents an anomalously big surface magnetization of order one, and its decay with L will be governed by the measure 
of these rare samples as a function of L. However, since the limit distributions above are a priori not expected to 
describe well the tail w ^ 1/VL outside the scaling region, and since the approximations (|38ll are not valid anymore 
for ~ 1, it is necessary to discuss separately the samples having Zl ^ 1 from the starting point 124|l of the saddle- 
point analysis. For these samples, the minimum (—El) of the potential U{x) on the interval [0, L] has to remain of 
order one. At the critical point /i = 0, the decay of the mean surface magnetization in the canonical ensemble pj, 

[ml]cano OC ^ (45) 

corresponds to the scaling of the probability for a Brownian path to remain above its starting point during a length 
L. In the microcanonical ensemble, the presence of the different decay 0,01 

[ml]micro OC ^ (46) 
L — ^oo ±j 

can be understood as the ratio between (i) the scaling as l/L^/^ for the probability of the first return to the origin 
after a length L, (ii) the scaling as for the probability to be at the origin at L. These two scalings, reflecting the 

measure of the rare samples having Zl ^ 1, can actually be recovered from the scaling functions H43|l as computed 
in 0, but the numerical prefactors given in equations (12-13) of Q should not be considered as exact. Actually 
these prefactors are not expected to be universal, as was discussed for other quantities such as the average end-to-end 
correlation in . 

In conclusion of this Section, we have shown that the limit distributions for the log of the surface magnetization 
obtained in ^ can be derived via a saddle point method in each sample, and have thus a very simple probabilistic 
interpretation : they describe the probability distribution of the minimum of a finite-size biased random walk, with 
different boundary conditions in the canonical and in the microcanonical ensembles. These limit distributions for 
the log of Kesten random variables should be considered as the analog of the Central-limit theorem for the log of 
a product of random variables ||33 : they describe well typical samples and are expected to be universal. However, 
the computation of some averaged quantities which are dominated by rare events cannot be computed exactly from 
these limit distributions in the scaling regime. In the two next Section, we thus compute exactly the average surface 
magnetization in the two ensembles. 

IV. EXACT RESULTS FOR THE SURFACE MAGNETIZATION IN THE CANONICAL ENSEMBLE 
A. Exact expression for the averaged surface magnetization 

The probability distribution of the random variable Zl H15|I in the case where the process {F(x)} is a Brownian 
motion has been already much studied by various methods I25ll3lll3^ . Here, the most convenient starting point 
is the following exact result for its Laplace transform derived in [2^ 



r(i + n)r(i + ^-n) 



[6 \cano — / ^ 6 



0<n<f 

f +C30 



H I dqe^^^^ ^gsinhTrg 

271-^ Jo 




(47) 
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where fj, is the exponent (|22|l . and where a represents the strength of the disorder (|21|l . We refer the interested 
reader to the previous studies H^l^ for a detailed discussion of this result. Here, we are interested into the surface 
magnetization mf which can be computed from the result (|47|l via the identity 



[ml]cano = [(1 + ZL)~^/%ano = ^ / 

Vi" Jo 

which yields in terms of the the Whittaker function (|B3p 



-j-oo 



[ iJcano ^ > Y{l + n)T{l + ^i-n) 



2 ' 2 



+ 00 



47r5/2 Jo 

As a comparison, the negative moment of order (1/2) reads 







2 




-1^ 




2 r 

















2' 2 1 cr 



(72 



0<n<-!i 



+- 



+00 



r(i + 7i)r(i + ^-n) 







2 





















47r5/2 7o 

i.e. the two results are very similar : the only differences are in the the factors inside [..] in (|49fl . 



(48) 



(49) 



(50) 



B. Ordered phase ^ = —25 > 

For /i > 0, there exists a limit distribution in the large L limit, whose generating function is simply given by the 
term n = in H47|l 



r(/i) 



It corresponds to the probability distribution [22, |2 



r(^) \gZ,, 



(51) 



(52) 



For the RTFIC, the case > indeed corresponds to the ordered phase 5 < Q where the surface magnetization 
remains finite in the thermodynamic limit 



i(oo)]cano(M > 0) = ^JL^/G ■ 



;r(i + A^)) 



r(i + M) 



(T 2 £2^ VV_li ii 

2'2 \ cr 



As a comparison, we have 



r(i + M) 



(53) 



(54) 



[Z oo' ]cano — 

C. Critical point = 

At the critical point /i = 0, the result 1471) for the Laplace transform of the probability distribution of Zj^ simplifies 
into 



[e-^^%ano = - / dqe-^'i" cosh^|if,J 2,/^ 



(55) 
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and the result (|49|l for the average surface magnetization becomes using (|B4I 



,{^^ = 0) = — dqe-^^^K,. ] (56) 



TT /n ^ \ 2a 



£2. /I 



^/naL \2cr 
whereas the negative moment of order (1/2) reads 



K'^ (^) 



0(1/l3 



(57) 



[Z^^Uno{f^ = 0) = ^ (58) 



Both results are indeed dominated by rare configurations of measure 1/\/L which have a of order 1, in agreement 
with the previous studies P, . The different prefactors in the leading terms of (|57I58|I can be understood via the 
series expansion 

(1 + Zl)-'^^ = Zl^l^ - {ll2)Zl'"^ + ... (59) 

Since the averages are dominated by the rare events where Zl is of order 1 with probability of order 1/\/L, all the 
terms in the series expansion will contribute to the prefactor found in l|57|) . 

It is actually interesting to consider the distribution of Z^ among these rare events : the Laplace transform has for 
leading term as i ^ oo with p fixed 

\e'^^-Uno,^ -^kJ2.[^+o(-^\ (60) 



which corresponds after Laplace inversion to the following behavior for the probability distribution -Pl(z) at large L 
with fixed z 

This result describes the tail of rare events with z fixed outside the scaling region Inz ~ VL studied in Section Hill 



D. Finite-size scaling function in the critical region 



In the critical region parameterized by the rescaled parameter 7 (|40|l . we obtain that the distribution of Zl of order 
one is characterized by the Laplace transform 



c can. 



6'(7 > 0)7^A + l 



I7l 



+ 00 



O 



(62) 



where the term containing the theta function comes from the bound state n = 0, and the other terms from the 
continuum. Actually, we may rewrite for arbitrary sign of 7 



\p-pZL 

|C- Jean 



0(7) 



VTnrL 



+ 7 



+00 



dve " 



O 



(63) 



This result is thus completely factorized into (i) the distribution of z among the rare events at the critical point H60(l 
times (ii) the function of 7 inside [..] that coincides with eq (12) of and that represents the probability of = in 
the scaling function discussed in previous Section IlIII As a consequence, the finite-size scaling form for the averaged 
magnetization reads 



^0(7) 



V TTOL 



e -1 



+00 



dve-'" 



0|i 



(64) 



where is a non- universal constant, found here to be c™ = e"i^ (^) (|57|) . and where all dependence in 7 
corresponds to the prediction of the scaling regime as derived in 0. We now compute the same observables for the 
microcanonical ensemble. 
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V. EXACT RESULTS FOR THE SURFACE MAGNETIZATION IN THE MICROCANONICAL 

ENSEMBLE 

In this Section, we derive exact results for the continuous Kesten variable with a Gaussian disorder in the presence 
of the niicrocanonical constraint 

/ dyF{y) ^ FoL ^ fiaL (65) 
Jo 

and we compare with the results of the canonical ensemble given in the previous Section. 

A. Probability distribution 

The Laplace transform of the probability distribution of Zl in the niicrocanonical ensemble may be written as 
follows in terms of path-integrals 

JVF{x)e-^loiFi-)-Fof-pIo ^""'"'j (/^^ dxF{x) - FoL) 

JVF{x)e-^lo[F(-)~FoVg(^jL^^p^^^^p^^'^ 

where the delta function represents the microcanonical constraint (ESJ. Rewriting this path- integral over the force 
F{x) as a path-integral over the potential U{x) — dyF{y) as in |25j, we obtain 

IuiQ)=o T^U{x)e 4. Jo I j 
The denominator is a simple Gaussian propagator 

VU{x)e~'^ ^0 [^n-) ^ P 4FI = (68) 

U{o)=o vAnaL yAnaL 

The path-integral in the numerator is equivalent to a propagator of quantum mechanics corresponding to the 
Hamiltonian H = —ad^/dU^ +pe^^ . Expanding this path-integral into the eigenstates ij^kiU) of the Hamiltonian as 
described in (2^, we finally get with Fq — 



-"^'^Ucro = %/4^e^'^^ / -^e-'='^^^(0)Vfe(M^i) 



2 



— 00 

+00 



27r 



- ^3/2 V aLe-""^ j dqe--''''qsmhnqK^q y^]J ^ j K^q 1^2^ ^e' — j (69) 
which should be compared with H47|l . 

B. Surface magnetization 

The surface magnetization mf which can then be computed via the identity 

[mlU^ro = [(1 + ZL)-'/^Ucro = ^ / dpp-'/^e-P[e-P^^U,ro (70) 

Jo 



which yields 



^V^e^'^'^ £ dqe^T-igsinh^q^^ dpp-^/^e-PR.g (^2^^ (2^6-^^ (71) 
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C. Ordered phase /i — —25 > 

For /i > 0, we have seen in the canonical ensemble that there exists a limit distribution in the large L limit, given 
by H51|l . To recover this result in the microcanonical ensemble, we may use (IBlll at lowest order to obtain 



2i sinh irq 



which yields 



I"^ imicro 



r(l-zg) 



2jP 



T{l + iq) 



(72) 



-o-L 



-ah 



T{l-iq) 



(73) 



where the final limit is obtained by taking into account the saddles in the complex plane aX q — respectively for 
the two terms. The limit distribution in the thermodynamic limit is thus the same as in the canonical ensemble (|51ll . 



D. Critical point fi = 



At the critical point /i = 0, the probability distribution of Zl is characterized by the Laplace transform 



\rmcr< 



o(m = o) = 



r3/2 



+ 00 



dqe 4 1 qsinhnqKfg — 



which should be compared with H55|l . The surface magnetization H71|) reads 

[mlUicro= dqe- "^"^q sinh irq dpp-^^^e-PKfJ2J- 

As a comparison, the negative moment of order (1/2) reads 

[ZJ^^^\nicro{^J^ = 0) = TraVi / dqe^'' 'r'q tanh(7rq) 

Jo 



l-^ + 0(l/L^) 



(74) 

(75) 
(76) 

(77) 



Both are indeed dominated by the rare events where Zl ~ 1, which have a measure of order 1/i in the microcanonical 
ensemble, in agreement with previous studies 

It is now interesting to consider the distribution of Z^ for these rare events : the Laplace transform has for leading 
term as L — > cx) with p fixed 



L^oa (tL \ \ (J J \^ 

This leading behavior is actually very simply related to the analog result (|6()|l in the canonical ensemble 



-pZ L^rare 



7r([e 



IcanoJ 



(78) 



(79) 



which corresponds after Laplace inversion to a convolution for the probability distribution P£''"°(z) at large L with 
fixed z (|^ 



P 'micro 
L 



{z) ^ ^Ko(-)e 

L^ao aLz \ az 



(80) 



This result describes the tail of rare events with z fixed in the microcanonical ensemble, outside the scaling region 
In z '--^ studied in Section IIIII 
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E. Finite size scaling function in the critical region 

In the critical region inside the ordered phase n — —26 > 0, parameterized by the rescaled parameter 7 H40|l . the 
analog of H72|l reads with q = k/VoL 



yVa-L 



2i sinhTrg 



m-iq) 



which yields for the generating function with fixed Z at large L 



-pZl 



,(7 = ^VctL > 0) = 7- 



r{l + iq) 



'cfL V V 



The finite-size scaling for the average surface magnetization is thus simply 

7Cm 



,(7 > 0) 



(81) 



(82) 



(83) 



that should be compared with the corresponding result in the canonical ensemble Hfi4|l : c™ is the same constant, and 
the scaling function of 7 has been simply replaced by the factor 7 alone, that again corresponds to the prediction of 
the scaling regime as derived in Here, the amplitude of the leading term in \/\/L vanishes at the critical point 
7 = 0, as it should to recover the scaling 1/L at the critical point. 

In the critical region inside the disordered phase /i = —25 < 0, parameterized by the rescaled parameter 7 H4U|) . the 
argument in the Bessel function is now exponentially large in -v/L instead of exponentially small as in (|81(l . Taking 
into account the asymptotic behavior of Bessel function at large argument ljB2|l we thus obtain that 7 is not the 
appropriate scaling variable near the critical point on the disordered side. Instead, to obtain the finite-size scaling 
function, we need the scaling variable 



that leads to 



p — jiuL 



,{p = fiaL<0) ~ ±Ko(2.^]K„ 



P -p 
' -e 2 



After Laplace inversion, the probability distribution Pl{z) at large L with fixed z and p is given by 



L 



1 



du 



L^oo aLz Jq u{1 — u) 



O 



(84) 



(85) 



(86) 



So here, in contrast with the corresponding result (|53|) in the canonical ensemble, this result is not factorized into a 
function of z times a function of p. The finite-size scaling form for the average surface magnetization is thus given by 



ro{p = fJ-crL < 0) 



+0 



(87) 



In conclusion, the critical regime on the disordered side is thus governed by the scaling variable p — pcrL < in 
the microcanonical ensemble, whereas it is governed by the scaling variable 7 — p\J aL < in the canonical ensemble. 
To understand these results, it is thus useful to recall the origin of the presence of two different correlation length 
exponents in the RTFIC. 

F. Discussion on the two correlation length exponents 

The presence of two different correlation length exponents i' = 2 and = 1 in the RTFIC has been discussed in 
detail in The definition of these two length scale for the Brownian process {U{L) — C/(0)) can be summarized as 
follows 7] : the first length scale corresponds to the length ^ where the mean value < U{L) — C/(0) >= FqL — paL 
is of order one, which yields 



with i) — 1 
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The second length scale corresponds to the length ^ where most of the samples indeed have {U{L) — U (0)) of the 
same sign of the mean value, i.e. the scale ^/oL of the fluctuations should be of the same order of the mean value, 
which yields 

^ L with i^^2 (89) 

These two length scales appear in various quantities in the RTFIC. In particular, the averaged correlation is governed 
by the exponent v = 2, whereas the typical correlation is governed by the exponent f' = 1 0. More generally, 
the presence of these two length-scales is also well known in related models, for instance in the eigenstates of the 
Fokker-Planck operator in a Sinai potential |22l and in one-dimensional random- hopping Hamiltonian for fcrmions 

In view of this knowledge, the interpretation of our results obtained in this Section for the averaged surface mag- 
netization is as follows : in the canonical ensemble, the critical behavior is governed by the exponent v — 2, whereas 
in the microcanonical ensemble, there are two different exponents in the two sides of the critical point : the ordered 
phase is governed by the exponent v = 2, with an accidental vanishing of the amplitude at the critical point, as 
stressed in , but the disordered phase is governed by the exponent v = 1, in agreement with the analysis given in 
and in contrast with the interpretation given in . More generally, the appearance of these two exponents for some 
given observable in the two ensembles should not be too surprising : indeed, the requirement that almost all samples 
of size L indeed "know" the sign of Fq involves the length-scale ^ (I89|l in the canonical ensemble, whereas it involves 
the length ^ (|88|l in the microcanonical ensemble. So the very definitions of the two length scales ^ and ^ show that 
the microcanonical constraint can indeed play an important role for the critical exponents. 



G. Obtaining the exponent v from the scaling regime 



As a final remark, it is now useful to discuss how the exponent i> ^ I, found in this Section by an exact path-integral 
method for the averaged surface magnetization, can be understood from the knowledge of the scaling regime H43|l . 
that a priori only contains the exponent v ^ 2. The evaluation of the averaged surface magnetization from the scaling 
function in the disordered phase 

K]rcr(M < 0) = / dwQr^^^ro {w; j) e"^- (90) 
Jo 

shows that the origin of the "anomalous" exponent i> = 1 can be traced back to the presence of the theta function 
d{w > —7) in the scaling function Qmicro {w, 7), i-e. the probability density Qmicro (w; 7) vanishes at a finite positive 
value in w in the disordered phase. This remark will be useful in the next Sections on more complicated observables, 
where we will compute the probability distributions in the scaling regime, and where direct methods to compute 
exactly the averaged values such as the path-integral method used in this Section will not be available. 

VI. CORRELATION BETWEEN THE TWO SURFACE MAGNETIZATIONS IN THE TWO 

ENSEMBLES 

In this Section, we consider the correlation {mfrnD © between the two surface magnetizations that has been 
numerically studied in Q. We compute its universal limit distributions at large L in the two ensembles. 

A. Saddle-point method in each sample 

The correlation (mfrn^) can be expressed in closed form in an arbitrary sample as the surface magnetization lO 

mlml = [1 + Zl]"'/' [1 + Zl] (91) 
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where Zj^ and Zj^ are two Kesten variables defined in terms of the same random potential Uix) = F{y)dy in the 
continuous version (|15|l 

Zl = /^da;e^(")-^(^) (92) 







Zl = /'^da;e^(-)-^(^) (93) 







The saddle-point method in a given sample 124|) on the associated random potential U{x) = dyF{y) leads to 

Zl ^ e^(°)-^— Z,,„e, (94) 
Zl e^— (95) 



L- 



in terms of the minimum Umin and maximum Umax reached by the process U{x) on the interval [0, L]. Since both 
Kesten variables are typically large in the critical region, we may replace as in (I38II 

ln(m>^) ~ --XuZl- -Zl ^ — — (96) 

SO the limit distribution of ln(mfm|^) can be computed in terms of the distribution of the random variable 

Dl = A{L) - U{L) (97) 
where A{L) is the amplitude of the Brownian trajectory for Q < x < L 

A(L) = Umax ~ Umin (98) 

and where C/(0) = is the starting point and U{L) the final point. 

B. Joint distribution of the amplitude A and the end-point U 

The joint probability P'^^ [A^ U) of the amplitude Al = Umax — Umin Ub — Ua and the end-point U when starting 
at J7o = has for Laplace transform with respect to L (see Appendix IX| 

r+oo 

P^^'HA,U]p)= dLe-P^ p[^\A,U) 



where 



"'0 

nt A ^ ,TT,\ jic/ COshg[/(g(A- |£/|)cothgA- 1) +sinhg|[/|cothgyl 

0(A>\U\)e^ r-ri — 1 (.^^) 

a smh qA 



(100) 



The Laplace inversion yields 



-2k{k - 1){\U\ - 2kA)e'''^] (101) 
We now discuss the microcanonical ensemble and the canonical ensemble respectively. 
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C. Microcanonical ensemble 



In the microcanonical ensemble with the constraint U (L) — ^aL, the probability distribution of the amplitude 
Al = U„iax ~ Ujnin ^ Ub - Ua rcads (llUlfl 



^ ,\ ' ^ — = e{A > \h\<tL)—— > e 



{A - \^i\aL)k^ ' ' 



- 2 e 



2k{k- l){\fj.\aL-2kA)e 



kA\fj.\ 



\ ctL 

Since D = A ~ fiaL, we may now write the final results for the various cases /i = 0, /i > and /i < 0. 



(102) 



1. Critical case /j, — 
In the critical case, the probability distribution of the variable D = A reads 

+ 00 / 7-,2 



p(^=°)(i^) ^9{D>0)—ye{2k^—-S 

rr r, ' ^ \ rrJ, 



ah 



k=\ 



We thus obtain that the distribution of the variable ln(mfm|^) (|96|l takes the scaling form 



[7'L(-ln(m>i))], 



' micro ' 



—2 ln(mf m|^) 
Vol 



with the scaling function 



+ 00 



k=l 

+ C30 



2^2 



2n^7r 



3 e-" ^ 



(103) 



(104) 



(105) 



(106) 



These two series representations are useful to study the asymptotic behaviors at small and large arguments respectively. 
The decay at large V is given by 



PniicroC^^ — 8X^ 6 



and at the origin, there is an essential singularity 



p (V) — 



(107) 



(108) 



This behavior of the scaling function at small argument yields the following dependence in L for the averaged corre- 
lation of the two surface magnetizations 



r s SI scalinq 



+ 00 



icro 

{V)e 



L — ^cjo 



^2 



(109) 



that should be compared with the numerical fits discussed in 



2. Ordered phase /i > 
In the ordered phase /i > 0, there exists a limit probability distribution for D a,s L oo 

Pt^°\D)\ =e{D> 0)fi^De~''° 



(110) 



This result corresponds to the convolution of two independent barriers against the bias as it should. 
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3. Disordered phase fj, < 

For ^ < 0, the minimal value of D is -Dmm ~ —2^aL, so it is convenient to set 

D = 2\fi\c7L + W 

and there exists a limit distribution for W in the large L limit that reads 



(m<o) 



(W) 



(111) 



(112) 



This results simply represents the distribution of the sum of two barriers against the bias : the comparison with the 
result H37|l concerning one surface magnetization shows that it can be interpreted in the same way. 

4- Finite-size scaling function 

Finally, in the vicinity of the critical point, we may write the following finite-size scaling form for the probability 
distribution of (— In mfm^) 



Pl{— IrimlmD 



--Q[V 



-2 In mfmr , 

^;7 = ^Vo-L 



'aL 



(113) 



-t-oo 



where the scaling function reads 

k— — OQ 

k^V + 7 - |7|) i2kV + {2k - e'=''(^+'^' + 2k{k - l){2k(V + 7) - |/i|cri)e'=l'^l(^+^) 



(114) 



By comparison with the results obtained before for the surface magnetization, we expect that the averaged corre- 
lation of the two surface magnetizations is governed by the exponent v — 2 in the ordered phase and by the exponent 
i> = 1 in the disordered phase, for the reasons explained around Equation (|90|l : the scaling function 1)1141) is defined 
in terms of the exponent 1/ = 2, but the presence of the theta function d{T> > —27) that forbids a finite interval in V 
near the origin in the disordered phase yields that the averaged value [mlm^]raicro will be governed by the exponent 
i> = 1. 

D. Canonical ensemble 

For the canonical ensemble, it is more convenient to characterize the probability distribution of D by its Laplace 
with respect to L using the result l|99|) 



p('''(Z?;p) 

= e{D > 0)D 



+00 



dLe'P^ 



cano jQ 

-\-oo 



dA / dUP^''\A,U;p)S{D-{A-U)) 



]j/^_i\qD coshqDv coshqD{v — 1) — smhqD 



11 (jsinh qDv 

+e{D>0)D t ^^'^^^ ~ ^) coshgj:>« coshgj:>(l - v) + sinh qD{\ - 2v) 

J 1/2 (Tsinh^qDu 



In particular at criticality fi — 0, the distribution of the variable D takes the scaling form 



1 .-pio) (v= ^ 

- ' catio I 



where the scaling function is determined by the integral transform 

dV 
p2 



dV _4„(o) /■+°° , q cosh qv cosh q{v - 1) - sinhq 



2 sinh qv 

^ , (7(21; — 1) cosh cosh 5(1 — w) — sinh (7(2^ — 1) 

dv ; — ^ 

1/2 2 sinh qv 



(115) 
(116) 

(117) 



(118) 
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In particular, there is a logarithmic divergence in the limit q ^ 

r§e-^Vi±iV) ^ llni (119) 
Jo l^'^ 9-^0 Z q 

that yields the linear behavior near the origin 

'P^loi'D)^c^V (120) 

in contrast with the essential singularity found above for the microcanonical ensemble H108|) . As a consequence, the 
averaged correlation of the two surface magnetizations decays algebraically in L 

r+co 1 

{^lrnl]Z:r^ dPP_„(P)e-^^ « j (121) 

in contrast with the stretched exponential behavior found for the microcanonical ensemble (I109|l . This big difference 
between the two ensembles may be understood as follows : in the microcanonical ensemble where U{L) = 0, the 
requirement to have a small D is equivalent to the requirement to have a small amplitude A — Umax ~ Umim which 
is a very strong constraint for a Brownian trajectory that leads to an essential singularity for the probability of small 
D. On the contrary, in the canonical ensemble, the requirement to have a small D = A — U{L) is not equivalent 
to the requirement of a small amplitude A : the freedom in U{L) allows to have a large amplitude A provided that 
U{L) is of the same order : the samples having a small D are thus those having Umin ~ ( probability 1/^/L) and 
Umax ^ U{L) ( probability \/\/L) and thus the measure of these samples is of order 1/i, in agreement with the above 
result lfnT|l . 

Finally, we note that the decay as 1/L for the averaged correlation of the two surface magnetizations (II 2 1|) in 
the canonical ensemble is similar to the decay found in [iSj for the averaged spin-spin end-to-end correlation in the 
canonical ensemble. We will discuss the relation between these two observables in more details in the next Section. 



VII. GAP AND END-TO-END SPIN-SPIN CORRELATION IN THE TWO ENSEMBLES 

In this Section, we consider the end-to-end spin-spin correlation and the gap in finite samples. In contrast with 
the surface magnetizations studied in previous Sections, these observables cannot be written in closed form in terms 
of all the random couplings. However, they have been studied in details in the canonical ensemble via the real-space 
renormalization approach 13], whose predictions are in good agreement with the numerical data [3. Il3|. The aim in 
this Section is to derive the corresponding results for the microcanonical ensemble. 



A. Gap and End-to-end spin-spin correlation in a given sample 

We refer the reader to the papers 01^3 for a detailed description of the real-space renormalization approach. Here, 
we will only recall its results concerning the interpretation of the end-to-end spin-spin correlation C{L) Q and of the 
gap A(L) © for a given sample jl3j| : for a given realization of the random potential U{x) for < x < L, the gap 
A(L) is asymptotically determined by the last renormalized ascending barrier G via 

G = -lnA(i) (122) 
and the correlation C{L) is asymptotically determined by the variable A = G — U (L) + U (0) via 

A = -lnC(i) (123) 

In particular, the RSRG predicts a very simple relation between the gap and the correlation in each given sample 

G-A = U{L)-U{0) (124) 

This prediction has been confirmed by the numerical studies at the critical point, in the canonical sample where 
(G — A) / VL was found to converge towards a Gaussian distribution |3, 0| , as well as in the microcanonical sample 
where (G — A)/-\/L was found to converge towards a delta distribution 0- 

We refer the reader to 0,^3 for a full theoretical justification of these results. However, it is interesting to mention 
here that the RSRG prediction for G{L) in the presence of a fixed boundary condition exactly corresponds to the 
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saddle-point analysis of Section ITTTI for the surface magnetizations. Indeed, if the spin ctl = 1 is fixed, it cannot be 
decimated, and the last ascending barrier is constrained to contain this end-spin : the variable G is thus determined 
by the minimum via G — U{L) — Umin and the variable A reads A = Umin ~ U{Q), in agreement with the saddle-point 
analysis of Section ITTll for the surface magnetization m\. When the boundary spins are both free, the gap and the 
correlation arc not in general simple functions of the end-values U{Q) and U{L) and of the extrema Umin and Umax 
of the random potential, but are determined by the properties of the last stage of the renormalization procedure as 
we now recall. 



B. Last stage of the renormalization procedure 

In the RSRG approach 0, one needs at a given scale F the measure B^{F, I) of the ascending (resp. descending) 
Brownian paths that goes from U{Q) to U{1) = F (resp. U{1) = —F) with F >r and with no return of more than F, 
i.e. for two arbitrary points < xi < a;2 < ^, the potential has to satisfy U (a^i) — U{x2) < F (resp. U{x2) — U{xi) < F) 
[35I l36l| . In Laplace transform with respect to I these measure read Q 



B^{F-p)= / dle~P'B^{F,l) =0(i^-r)e±^r^-^e-(^-i^)(T^+«™"^«r) (125) 
JO sinh(/l 



where q = yjp + and cr = 1 to simplify the notations in this Section. For the boundaries, one also needs the 
measure E^{F,l) of the ascending (resp. descending) Brownian paths that goes from C/(0) to U{1) = F where F >0 
is the maximum (resp. U(l) = —F is the minimum) with no return of more than F. In Laplace transform with respect 
to I these measure read 0, I^E 

Ep{F;p) ^6'(i^)e-^(=Ff-H9cothgr) ^^26) 

The joint distribution Pl{G, A,U) of the variables G = — \nA(L), A = — In C(L) and of the end-point Ul = U 
(with U (0) = 0) is then determined by the last stage of the RSRG as 

PUG, A, [/) = J dhdl2dhdARdALE^{KR)B+{G)E^ {Al)5{L - [h + I2 + h))SiA - {Al + Ar))S{U - (G - A)!jl27) 

i.e. in Laplace with respect to L 

P(G, A, U;p) ^ e{G)9{A)Ae-^'^^+^'°'^''^^e^^ . ^ MU - [G - A)) (128) 

smh qG 



C. Canonical ensemble 



In the canonical ensemble, the joint distribution P^'^\g, A) has for Laplace transform 



+00 



dLe^P^ 



Pi^^ (G, A) = / dUPiG, A, U; p) = e{G)0{A) 



qA 



sinh qG 



(G-A)g-A9coth9G 



(129) 



which indeed coincides with the result (45) of [13. We refer to [I^l for a detailed analysis of this result, and only 
quote in the following some important results in order to compare them with the microcanonical case below. 



1. At criticality fi = 
At criticality, the law for the gap alone has for Laplace transform 



+00 



dLe-P^ 



Pr\G) ^e{G) :""'V" (130) 



sinh Y^G 



y/pcosh y/pG 

that corresponds after Laplace inversion to the series given in Eq (57) in [3. In particular, the probability distribution 
of the rescaled variable g = G/\/X presents the essential singularity at the origin 

^cT„o(ff)^«„^e-^ (131) 
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that determines the decay with L of the averaged gap |l3| 

i^WTcanr' « Li/6e"K^^)'^' (132) 
On the other hand, the law for the correlation alone has for Laplace transform at criticality 

/' + OC 



4-oc 

dLe-P^ 







p|^=°)(A) =0(A)v^a/ dG^— — -e-^v1^™*'^^^ = 0(A)Aifo(VpA) (133) 

J cano jQ Smh y/pCr 

which corresponds after Laplace inversion to the simple result for the rescaled variable A = A/\/L 

Kal'oW^ei\)^e''^ (134) 
in excellent with numerical data [T^. As a consequence, the averaged correlation decays as a power of L 

L — ►oo L/ 

So in the canonical ensemble, at criticality, the gap and the spin-spin correlation have very different behaviors. In 
particular there is an essential singularity at the origin only for the gap and not for the correlation, and thus the 
decays for the averaged values are very different. 

2. Ordered phase fj. = —2S > 

In the ordered phase, the variable A is simply the sum of two independent barriers (U{0) — Umin) and {Umax ~ U{L)) 
against the bias [Ijl 

Pieo°^(A)c.^2Ae-''^ (136) 
whereas the variable G = U{L) + A is asymptotically distributed with the Gaussian distribution of U{L) alone 

^iL>°HG) ^ ^e-^ (137) 

3. Disordered phase /i = — 2<5 < 

In the disordered phase, the probability distribution of the variable G reads |0 

PiL<°)(G) (V)|He-IHG,-(L.^)eH.io ^^33^ 

that can be interpreted [T^ as the distribution of the maximal barrier among N ^ Ly? independent variables drawn 
with the exponential distribution |/i|e"l'^l'-^ of barriers against the drift, whereas the variable A = —U{L) — G is 
asymptotically distributed with the Gaussian distribution of U{L) alone jl^ 

Pi,^<°)(A) ^ -_=e-Mi^ (139) 

D. Microcanonical ensemble 

1. At criticality fi — Q 

At criticality /i = 0, we obtain that the the joint distribution pj^'^~^\G, A) is determined by the Laplace transform 



+00 







p1^=°)(g,a) 



P(G, A, U = 0;p) = 5{A - G) .-^^""'^^^ (140) 

smh Y^G 
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So here G and A are identical, and the rescaled variable 



9^^ = ^ (141) 



is distributed with the distribution Vmicro{g) determined by the integral transform 



^ / ^V^.cMe-f' = _^e-^-*''^ (142) 
Vtt Jo 9 smh ^ 



In particular, the following divergence near q 



Vmicroigy ^ OC ; 5-6^+^ (143) 

yields the following essential singularity for the probability distribution Vmicroig) at small argument 

Vrmcroig) OC ^6^^ + " (144) 

This behavior at small 9 of the scaling function characterizes the rare events having a gap Al = e~^^ (or equivalently 
a correlation Cl — e~^^) of order one, which are expected to dominate the mean values of the gap and of the 
correlation. The dependence in L of these mean values may thus be estimated by a saddle point method which yields 

/.+00 

■\scaling _ iscalmg / J„-n f^^^-QvL 



[^lZZ7 = [C^l]_:7 - dgVrn.cro{9)e-^^^ (145) 

0. £-l/6e-i(2.^L)^'^+2(2.^L)Ve ^^^^^ 
L — >oo 

So here we obtain that there is a sub-leading power L^^^ in the exponential with respect to the leading decay as 
E~'^^ in contrast with the canonical ensemble H132(l . 



2. Ordered phase p = — 25 > 



In the ordered phase, the variable A is again the sum of two independent barriers {U{0) — Umin) and {Umax — U{L)) 
against the bias, as in the canonical ensemble H136|l . This distribution of A now completely determines the probability 
distribution of the variable G = fiL + A, in contrast with the Gaussian distribution of the canonical ensemble H136(l . In 
particular, the inequality G > fiL indicates that the average of the gap A{L) = e~'^ will be governed by the exponent 
V = \ instead of the exponent v — 2, for the reasons given around Equation H90|) concerning the case of the averaged 
surface magnetization. 



3. Disordered phase /i = — 2<5 < 



In the disordered phase, the distribution of the variable G will be asymptotically the same as in the canonical 
ensemble H138|l . since the interpretation 13] as the the maximal barrier among N independent barriers against the 
drift still holds in the microcanonical ensemble. However, the distribution of the variable A is now completely 
determined by the relation A = + G, in contrast with the Gaussian distribution (|138|l of the canonical case. In 
particular, the inequality A > indicates that the average of the spin-spin correlation G{L) = e^^ will be governed 
by the exponent v = 1 instead of the exponent i/ = 2, for the reasons given around Equation H9U|) concerning the case 
of the averaged surface magnetization. 



E. Comparison between C(L) and mfmf^ at the critical point 

In the ordered phase, where the two end spins have a spontaneous magnetization, it is clear that the limit C(oo) 
corresponds to the product of two independent end-point magnetizations of an infinite system p^ . Moreover, it was 
found numerically 2] that the observables C{L) and {mlml) were actually still closely related at the critical point 
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: the limit distributions were found to be very close in the microcanonical ensemble, but different at large argument 
in the canonical ensemble 2J. On the other hand, the difference at small argument is much smaller in the canonical 
ensemble, and the average values [C{L)\cano and [m\m\\cano were found to be almost indistinguishable numerically 
in the canonical ensemble 2J. It is thus interesting to discuss these questions here. 

From the RSRG approach, it is clear that if the positions Xmin and Xmax of the minimum Umin and Umax satisfy the 
order < Xmin < Xmax < L, then the largest ascending barrier is simply given by the amplitude G = Umax — Umin- 
Then the variable A for the spin-spin correlation is given by A = {Umax — U{L)) + {U{0) — Umin) and thus exactly 
coincides with the variable D = A{L) — U{L) ()97|l determining the correlation of the two surface magnetizations. 
In particular, in the ordered phase, all samples satisfy < Xmin < Xmax < L in the asymptotic limit L ~^ co, and 
the variables A — D are the sum of two independent barriers against the bias. At the critical point, half of the 
samples satisfy the constraint Xmin < ^maxj and thus we obtain that for one half of the samples at criticality, both 
in the canonical and in the microcanonical ensemble we have the identity A = D that yields the identity of the two 
rescaled variables (lnC(L))/-\/L = (lnmfm£)/-\/i- On the other hand, for the other half of the samples which have 
Xmax < Xmin, the largest ascending barrier is smaller than the amplitude G < Umax — Umin and thus A < D are 
two different variables. In conclusion, the identity A = D for half of the samples could explain the similarities found 
numerically in j2j for their probability distributions. 

Concerning the mean values, the interpretation given after Equation (|121|l for the decay as 1/L of the averaged 
correlation of the surface magnetizations in the canonical ensemble yields that the relevant rare events precisely satisfy 
the order < Xmin < Xmax ^ L. As a consequence, we expect that the leading term in 1/L of the averaged spin-spin 
correlation exactly coincides with the leading term in 1 /L of the averaged correlation of the surface magnetizations in 
the canonical ensemble, in agreement with the numerical finding that both quantities were almost indistinguishable 

Q 

VIII. CONCLUSIONS 

The example of the RTFIC considered in this paper shows that the question of the ensemble dependence for the 
critical properties of disordered systems is rich and instructive. Indeed, even if the two ensembles are expected to be 
equivalent in the thermodynamic limit, their finite-size properties may nevertheless be quite different. 

At criticality, even if the appropriate rescaled variables for all observables considered here are the same in both 
ensembles, the probability distributions may be quite different in the two ensembles. In particular, some asymptotic 
behavior may be given by two different power-laws in the two ensembles (as in the case of the surface magnetization), 
or even by a power law in one ensemble and by an essential singularity in the other ensemble (as in the case of the 
end-to-end correlations). As a consequence, the decay with L of averaged observables can be drastically different in 
the two ensembles : for instance, the end-to-end correlations decay as a power law in the canonical ensemble, but as 
a stretched exponential in the microcanonical ensemble. 

Off criticality, the probability distributions of rescaled variables involve the exponent i/ = 2 in both ensembles, but 
averaged observables are again governed by rare samples, whose measures can be very sensitive to the microcanonical 
constraint. As a consequence, the critical properties of a given averaged observable may be governed by two different 
exponents in the two ensembles. For instance, we have shown via an exact path-integral method that the averaged 
surface magnetization in the disordered phase was governed by the exponent v — 2 in the canonical ensemble and by 
the exponent v = 1 in the microcanonical ensemble, in agreement with the interpretation of 0, and in contrast with 
the interpretation of '5| . We have moreover concluded from the domains of definition of probability distributions of 
other observables in the scaling regime, that similarly, (a) in the disordered phase, the averaged correlation of the two 
surface magnetizations and the end-to-end spin-spin correlation are governed by the exponent = 2 in the canonical 
ensemble and by the exponent = 1 in the microcanonical ensemble (b) in the ordered phase, the averaged gap is 
governed by the exponent v ^ 2 in the canonical ensemble and by the exponent v = 1 in the microcanonical ensemble. 

As a final remark, let us mention the very recent preprint |37| , where the RSRG method is generalized to study the 
end-to-end energy-energy correlation Ce{L) (in the canonical ensemble), and where averaged observables are again 
governed by rare events. 
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APPENDIX A: JOINT DISTRIBUTION OF THE AMPLITUDE A AND THE END-POINT U 

The probability G[u^_ij^^{U, Uq; I) for a random walk with bias _fo = fia to go from Uq = to U during L in the 
presence of two absorbing boundaries at U — Ua < and U ^ Ui, > has for Laplace transform [35l l3^ 



nr ^\TT ^ _ ^Ji (u -u„) sinh q(U -Ua)smhq{Ub- Up) 
[Ua,Ut]^ 'P' °> ~ 17777177777 777 



^lUa.,U,] 



aqsiBhq{Ub - C/q) 



if Ua<U<Uo 



(UMUo) = et(^-^o) Binhg(t/.- L^) sinh,(L/o-L/.) < ^ < 

aqsmhq{Ub - Ua) 



where 



CT 4 



(Al) 
(A2) 

(A3) 



The joint probability of the end-point U , of the minimum value Ua and of the maximal value Ub when starting at 
f/o = has thus for Laplace transform 



P^'^Hu,Ua,Ub-,p) = -du^du,G\^%^^{U,p\Q) 

iLjj q sinhg[/sinh(7([/f, + Ua) — 2cosh(7J7sinh(7C7a smhqUb- 



— P 2 ^ _ 



sinh^ q{Ub - Ua) 



(A4) 
(A5) 



Finally, the joint probability of the amplitude Al — Umax — Umin = Ub — Ua and the end-point U when starting at 
Uq — has for Laplace transform 



P'^''\A,U;q) dUa I dUbe{Ua<0<Ub)e{Ua<U <Ub)P'^^{U,Ua,Ub)5{A-{Ub-Ua)) (A6) 



+ 00 



mm(0,C/) 



0{A > \U\) I dUaP'^^'HU, Ua, Ua + A) 

' max{ — A,U — A) 



This finally leads to the result H99|) given in the text. 



APPENDIX B: USEFUL PROPERTIES OF BESSEL FUNCTIONS 



(A7) 



The series expansion of Bessel function at small argument reads 



■E 



2i sinliTTO ^-^ kl 

^ k=0 



tq 



T{-iq + k + l) Tiiq + k + 1) 
whereas the asymptotic behavior at large argument is given by 

1 r{iy + k+^) 



Efc!(2z)fcr(i.-fc+i) 



A useful integral is 



(T i^ffJ.+ l + iq\^f^J. + l-iq 



JO \ \l <^ J 2 V 2 J \ 2 

where W is the Whittaker function. For the special case n — 0, it simplifies into 



2 ■ 2 
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